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ABSTRACT 


A South China Sea wann-core eddy with sea surface temperatures (SST) higher than 29.5 ° C, 
recently reported by Chu and Chang (1995), appears in the central South China Sea (west of Luzon 
Island) in boreal spring, and strengthens until the onset of the summer monsoon (mid-May), then 
weakens and disappears at the end of May. Although its size and intensity varies, the warm-core 
eddy releases large moisture and heat fluxes into the atmosphere and in turn affects the monsoon 
circulation. The transient features and interannual variabilities of the warm-core eddy have not yet 
been studied. In this study we use 189,059 temperature profiles from the U.S. Navy’s Master 
Oceanographic Observation Data Set (MOODS) to investigate the eddy features. We use Optimal 
Interpolation to create a ten-day interval synoptic data set for 1964-84 on a 0.5°x 1° grid (finer 
resolution in zonal direction) from the MOODS SST data. Then, we formulate an Emprical 
Orthogonal Functions (EOF) analysis. The first EOF mode accounts for 35.5% of the variance and 
is characterized by a field which is similar to the mean field. The second EOF mode accounts for 
21.4% of the variance and is characterized by a meso-scale eddy (116°-119°E, 16°-19°N) to the 
west of Luzon Island. The size of the eddy is around 100,000 kml The time series of EOF2 behave 
differently between El Nino and non-El Nino years; a warm-core eddy often appears during the non- 
El Nino years, and a cold-core eddy usually exists during the El Nino years. A new scenario of the 
air-sea feedback is proposed for the warm-core eddy formation in non-El Nino years and the cold- 
core eddy formation in El Nino years. A numerical simulation from a three dimensional primitive 
equation model shows the formation of a warm-core eddy under climatological forcing (more likely 
for the non-El Nino condition). 
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I. INTRODUCTION 


The South China Sea (SCS) has a bottom topography (Figure 1) that makes it a 
unique semi-enclosed ocean basin that is overlaid by a pronounced monsoon surface 
wind. Extensive continental shelves (less than 100 m deep) are found on the western and 
southern parts, while steep slopes with almost no shelves are found in the eastern part of 
SCS. The deepest water is confined to a bowl-type trench. The maximum depth is around 
5,000 m. 

Based on limited data sets, studies show that eddies in the SCS are predominantly 
cyclonic in winter and spring and anti-cyclonic in summer, with sizes ranging from small 
to meso-scale (Huang et al., 1994). Both cold and warm eddies exist in the SCS. Dale 
(1956) and Uda and Nakao (1974) reported a cold eddy off the central Vietnamese Coast 
in summer. Nitan (1970) found a cold eddy located to the northwest of Luzon. Reports 
from the South China Sea Institute of Oceanology (SCSIO, 1985) indicate that in the 
central SCS, a warm-core eddy appears in summer and winter, but more closer to 
Vietnam in summer at the surface. 

During boreal spring, the surface atmospheric circulation pattern changes from the 
pattern of the northern winter monsoon to that of the northern summer monsoon. In 
winter, a semi-permanent surface subtropical high pressure center prevails over the 
northwestern Pacific with its east-west axis located north of 20°N (Figure 2a). Cold 
surges tend to move directly from southeastern China to the SCS with north-northeasterly 
winds. The winter cold surges tend to enhance the atmospheric convective activity in the 
Malaysia-South China Sea Region (Cheang, 1980). 

During spring, the surface atmospheric subtropical high pressure center moves 
closer to the equatorial region, with its east-west axis between 15°N and 20 °N over the 
northern region of the SCS and northern Indochina Peninsula (Figure 2b). Cold surges do 
not move directly out from southeastern China. They travel a long distance across South 
Korea and southern Japan and around the subtropical high in the western Pacific. They 



arrive in the southern region of the SCS as northeasterlies. Such a semi-permanent 
atmospheric surface anticyclone results in downwelling in the upper ocean. 

Recently Chu and Chang (1995) reported a warm pool (SCSWP) in the central 
SCS (west of Luzon Island) during the late spring season, that strengthens until the onset 
of the summer monsoon (mid-May), then weakens and disappears at the end of May. 
Although its size and intensity varies, the SCSWP releases large moisture and heat fluxes 
into the atmosphere and in turn affects the monsoon circulation. 

The thermal feature of eddies (warm-core or cold-core) greatly impacts on the 
monsoon circulation. Therefore, it is very important to know what kind of eddies occupy 
the central SCS in winter and spring seasons. Can we identify these eddies from the most 
complete data set - the U.S. Navy’s Master Observational Oceanographic Data Set 
(MOODS)? What are the major features of the eddies? 

The data set that we have used is a subset of MOODS and includes 189,059 
profiles birmed for the desired area coverage. The data set covers the years 1929 through 
1994 and comes from numerous instrument types including XBT’s, CTD’s and reversing 
thermometers. Our study field covers the area 5°S - 25°N and 105°E - 120°N. The data 
were binned for every ten days creating 756 consecutive temporal groups called “points”. 
The Optimum Interpolation (01) technology was applied to create a gridded data set and 
the GMT Software (Wessel and Smith, 1993) was used to plot SST contours. 

After the gridded data set has been established, we divided our study into three 
parts. First, a composite analysis of the gridded SST data was done to obtain the 
ensemble mean SST field and the mean seasonal variation. The latter shows strong 
evidence of the existence of a warm pool in the central SCS before the summer monsoon 
onset (April and early May). Second, the transient features of the SST field were 
quantitatively investigated using the method of Emprical Orthogonal Function (EOF) 
analysis. An air-sea feedback hypothesis was proposed after the EOF analysis. Third, in 
order to investigate the mechanism of the SCS warm-core eddy formation during boreal 
winter and summer under climatological forcing and its possible effects on the summer 



monsoon onset, a primitive equation model with second-order turbulence closure 
(Blumberg and Meller, 1987) was used. 

This thesis is organized as follows: Chapter II depicts the data sources, including 
the U.S. Navy’s Master Oceanographic Observation Data Set (MOODS) and its temporal 
and spatial distributions as well as the establishment of gridded data set. Chapter III 
describes seasonal and interannual variabilities. Chapter IV discusses numerical 
simulations. 
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Figure 1. Geography and Isobaths in Meters, Showing the Bottom Topography of 
the South China Sea. 
















II. DATA SOURCES 


A. MASTER OCEANOGRAPfflC OBSERVATION DATA SET (MOODS) 

The MOODS is a compilation of observed ocean data worldwide consisting of (a) 
temperature-only profiles; (b) both temperature and salinity profiles, (c) sound-speed 
profiles, and (d) surface temperature (drifting buoy). These measurements are, in general, 
irregular in time and space. In this study, we analyze temperature profiles obtained during 
1964-84 fi-om a variety of instruments. Due to the shear size (more than six million 
profiles) and enormous influx of data to NAVOCEANO from various sources, quality 
control is a difficult task. Our study domain includes the area S^S to 25°N and 105° to 
120°E; the data set within this region consisted of 189,059 temperature profiles after 
rejecting certain data during quality control. The primary editing procedures included 
removal of profiles with obviously erroneous location, profiles with large spikes, and 
profiles displaying features that do not match the characteristics of surrounding profiles. 
In shallow water, this procedure can be partially automated but also involves subjective 
interpretation because of the undersampling of MOODS compared to the spatial and 
temporal variability of the oceanography. 

The temporal and spatial distribution of the MOODS data are irregular. The 
sparsity of profiles near the Sunda Shelf and western coastal region of Borneo is shown in 
Figure 3 and the temporal distribution from 1964 to 1984 is shown in Figure 4. Figure 4 
indicates that most observations were obtained during the Vietnamese War (1965-69). 
There are nearly 14,800 profiles (maximum number of observations in a year) in 1968 
(Figure 5). The minimum number of observations is in 1984 (near 3,000 profiles) (Figure 
6). These figures indicate that the number of observations within a given year exhibits a 
seasonal variation. 
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Figure 3. Spatial Distribution of MOODS Data during 1964-1984 
























Figure 4. Temporal Distribution of MOODS Data during 1964-1984 
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Figure 5. Temporal Distribution of MOODS Data for 1968 













































Figure 6. Temporal Distribution of MOODS Data for 1984 


The main limitation of the MOODS data is its irregular distribution in time and 
space. Certain periods and areas are over sampled while others lack enough observations 
to gain any meaningful insight. Vertical resolution and data quality are also highly 
variable depending much on instrument type and sampling expertise. There are two data 
sparse area: the Sunda Shelf area and the western coastal region of Borneo (Figure 3). 
There are frequent ten to 20 day gaps of no observations in the whole South China Sea. 
Spatial and temporal irregularities along with the lack of data in certain regions must be 
carefully weighted in order to avoid statistically-induced variability. 

B. ESTABLISHMENT OF GRIDDED DATA 

Our near ten-day interval data set is distributed evenly in time and imevenly in 
space. We use the Optimum Interpolation scheme to obtain a gridded data set. We use the 
i°x 1 ° monthly SST climatology of Levitus (1982) as a mean field. For the same time 
period, the SST perturbations V^ ^ for SST observations were obtained by 

subtracting the climatological value at the closest grid point from the observation. 
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Optimum interpolation assigns a weight to each observation that accounts for variation in 
spatial sampling. The interpolated temperature anomaly at the grid point, T'c, is a linear 
combination of the observed anomalies, 7’, ,...,7^^ with weights . 

Ta = 'ta,T, ( 1 ) 

The values of the weights, a, (; = 1,2,..., A^), are found by minimizing, in a least square 
sense, the difference between the interpolated value at the grid node and the true value 
there, and obtained from solving a set of N algebraic equations (Gardin, 1963) 

= (2) 

where i = l,2,...,N. Here A is the signal-to-noise ratio, is the autocorrelation 
between locations i and j , and is the autocorrelation between the grid node and 
location i. Here, the autocorrelation function (ACF) is defined by 

M = yir(r,)r(r,+r)dr, (3) 

where denotes the independent space vector defining the location of points in a 
sampling space R, and r is the space lag and the variance, ju is computed by paring 
the anomalies into bins depending upon their separation in space, r. The values of /u are 
obtained from calculating the correlation coefficient for all the anomaly pairs in each bin. 
The ACF will be constructed for the combination of different lags and in the isotropic 
case 


Tjim) = Tj{0)c}ip[-A^(mAry] 


( 4 ) 



where 77(772) denotes the ACF value in the bin with the horizontal separation 772A7-. At- = 
10 km is the increment for the space separation. A"' is the horizontal decorrelation 
scales. Therefore, the autocorrelations fj.-^ and are given by 

=ri{m,j), Hq.= 77(777^.) (5) 

where m.- and tttq, represent spatial separations between locations i and 7, and between 
location i and the grid node. 

The measured variance of the thermal fields is separated into signal and noise 

whereby 

5^ = s] + s‘ (6) 

The noise variance is brought on from two sources, geophysical and instrumentation 
errors. The geophysical error is unresolved thermal variability with scales smaller than 
the typical time and space scales between two temperature profiles, and the 
instrumentation error is negligible compared to the geophysical error (White and 
Bernstein, 1979). In calculating the ACF value of the first bin (zero to 10 kra lag) the 
temperature anomalies were never paired with themselves. The ACF value at the zero lag 

77(0) was extrapolated from the first bin into the origin, and therefore does not equal 1. 

The signal-to-noise ratio is computed by (Sprintall and Meyers, 1991). 



The larger the A, the less geophysical error exists. If 77(0) = 1, there is no noise, A = 00; 
and if 77(0) = 0, there is no signal, then A = 0. If A > 2, the ratio of the signal variance, 
sj, to the noise variance, s^, is greater than 4, which was considered a quite good 
condition by White et al. (1982) and Sprintall and Meyers (1991). 

For each grid node location, (1) results in a set of N linear algebraic equations to 
be solved for A unknowns, a,, by matrix inversion. In (1) the autocorrelations 
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and the signal-to-noise ratio X have been computed at the surface from the MOODS data 
set. Therefore, the N weights a, can be determined by solving (1), and thus any new 
observations , can be interpolated into any grid node G. 

We applied Optimal Interpolation to build up a 10-day interval synoptic data set 
for 1964-84 on a 0.5°x 1° grids (finer resolution in zonal direction) from the MOODS 
temperature data. Our gridded data set consists of 756 consecutive SST fields from 1-10 
January 1964 to 21-31 December 1984. An example of SST contours in 1980 is shown 
in Figure 7. 
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Figure 7. Consecutive SST Contours after the Establishment of Gridded Data in 
1980: (a) April 1-10, (b) April 11-20, (c) April 21-30, (d) May 1-10, (e) May 11-20, (f) 
May 21-31. The Solid (Dash) Lines Means SST Higher (Lower) than 28°C. The Dots 
Indicate the MOODS Stations. 
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III. SEASONAL AND INTERANNUAL VARIABILITIES 


A. SCS WARM-CORE EDDY CASE STUDIES 

Chu and Chang (1995) reported a warm-core eddy (SST >29°C) appears in the 
central SCS before the summer monsoon onset from analyzing the MOODS data for 
1966. A warm-core eddy often appears in spring and disappears in later May. The SST 
field in 1966 shows the following features: 

1. Flat SST Pattern in March 

This period is right after the cessation of the northern winter monsoon but some 
hint of the winter monsoon characteristics are still present (Figure 8a-c). Surface 
isotherms in March show a pattern of northeast-southwest oriented isotherms with a 
positive temperature gradient towards the southeast near the equator. The warmest water 
(>29°C) appears off the northwest coast of Borneo. SST decreases from 29°C near the 
Borneo coast to 24° at the southeast coast of China. 

2. Early Stage of Warm-Core Eddy 

This period is in the center of the transition season (Figure 8d). In early April, a 
warm pool with the size of near 500,000 km^ appears near the west coast of Luzon. The 
temperature of the warm water increases from 28°C (Figure 8d and e) to 30°C (Figure 
8f). 

3. Fully Developed Stage of SCS Warm-Core Eddy 

This period occurs in late spring before the summer monsoon in the SCS begins 
(Figure 8g). The warm water pool expands southwestward. During May 1-10, the central 
SCS was dominated by warm water (>30°C). A warm pool with SST higher than 30.5 °C 
is located at 111 °-l 16°E, 7°-12°N. This warm pool provides a surface heating source for 
the formation of the atmospheric low pressure center, which will contribute to the 
summer (southwest) monsoon onset. 
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Figure 8. South China Sea SST Evolution in 1966: (a) March 1-10, (b) March 11-20, 
(c) March 21-31, (d) April 1-10 (From Chu and Chang, 1995). 
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Figure 8. South China Sea SST Evolution in 1966: (e) April 11-20, (f) April 21-30, 
(g) May 1-10, (h) May 11-20 (From Chu and Chang, 1995). 
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4. Disappearing Stage of SCS Warm-Core Eddy 

This period is just after the summer monsoon onset (Figure 8h). The SST for May 
11-20 is relaxed back to the ensemble mean pattern and the SCSWP disappears. Warm 
water (>31°C) appears off the northwest coast of Borneo and decreases from Borneo to 
the southeast coast of China. 

5. Cold-Core Eddy Building-up by Coastal Upwelling 

Under the steady summer (southwest) monsoon forcing, cold water will be 
upwelled near the southeast coast of China and Vietnam. A cold-core eddy appears off 
the southeast China coast (<25°C) and off the Vietnamese coast (<28°C). 

B. COMPOSITE ANALYSIS 

The existence of the central SCS warm-core eddy is likely to be important for the 
air-sea heat and moisture exchanges and its effects on the summer monsoon onset. The 
warm-core eddy was reported to the west of Luzon Island during the winter-to-summer 
monsoon transition in 1966 (Chu and Chang, 1995). We now examine the data to see if 
we can obtain a warm-core eddy signal from historical data. To answer this question, we 
need to perform a composite analysis. After the gridded data is established, the SST data 
are represented by 7’(X(,y^,rj,f,), where (X(,y^.)are the horizontal grids, 1964, 
1965,... 1984, is the time sequence in years, and f,=l, 2, ... 36, the time sequence in ten 
day intervals. Before investigating the annual variation of SST, we define the following 
two temporal averages: 

= Ar = 1984 -1964 = 21 years (8) 

Atj 

which is the annual mean value, and 

nx>,yj)^^yT{x,,y^,t,), (9) 

which is the ensemble mean value. The ensemble mean (1964-1984) SST field over the 
SCS (Figure 9) shows the pattern of northeast-southwest oriented isotherms with a 



positive temperature gradient towards the southeast near the equator. The ensemble mean 
has a rather weak horizontal temperature gradient. T decreases from 29 °C near the 
Borneo coast to 26°C near the southeast China coast. 

The annual mean value relative to the ensemble mean, 
T{x-,yj,ti)= T {Xi ,yj ,ti) - T{x- ,yj) , shows the composite features of the annual mean 
SST anomalies (Figure 10). The annual mean SST anomaly (7) has the following 
features; 

1. Cooling Phase 

In March, T is negative almost everywhere throughout the whole SCS, and 
decreases from the west coast of Borneo to the southeast coast of China (Figure lOa-c). 
This indicates that the SCS is cooling during the winter and early spring seasons and that 
the rate of cooling is weakest near the Borneo coast and strongest near the Taiwan Strait 
(<-3“C). 

2. Early Warming Phase 

In early April (Figure 1 Od), positive T were found to the west of Luzon Island 
and the west coast of Borneo. At the Luzon Strait T is always negative until late April 
(21-30 April). This means that the early SCS wanning starts in the central SCS by some 
local effects (e.g., a warm eddy but not by the Kuroshio intrusion!) and at the south SCS. 
It also implies that the Kuroshio water entering the SCS in April might be cooler than the 
SCS water (at least at the surface). 

3. Central SCS Warming Phase 

In May (Figure lOg-h), the central SCS warm-core eddy becomes stronger than 
the early warming phase. During 11-20 May (Figure lOh), the central SCS warm-core 
eddy ( T > 1.8 °C) accounts for the maximum annual mean warm anomaly. 
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4. Mature Warming Phase 

After the summer monsoon onset, the June T becomes positive almost everj'where 
throughout the whole SCS, and increases from the west coast of Borneo («0.6°C) to the 
southeast coast of China («2.4“C). This indicates that after the summer monsoon onset, 
the SCS warming is weakest near the Borneo coast and strongest near the Taiwan Strait 
(>2.4°C). 

C. EMPRICAL ORTHOGONAL FUNCTION (EOF) ANALYSIS 

The composite analysis demonstrated the existence of a central SCS warm-core 
eddy in the aimual mean SST anomaly data. What are the transient features of this meso- 
scale eddy? What are the synoptic and interannual variabilities of the eddy? We can use 
the Emprical Orthogonal Function (EOF) technique to answer these questions. 

1. EOFs for the SST Anomalies 
The SST synoptic anomalies obtained by 

t{x„y^,r„t,)=T{x„yj,r„t,)-fix.„y^,t,) (10) 

are re-arranged into a NxP matrix, 7'(/-„,7 ^),m= 1,2,...,A; and p=l,2,...,P. Here P = 
756 is the total number of time points used for computing the covariance matrix, i.e., 21 
years of every 10-day measurements; N = 255, corresponds to the number of grids 
(i = 1,2,...,17;/ = 1,2,...15). Lorenz (1956) used the eigenvectors of the covariance matrix 
associated with a sample data field to calculate the time series of coefficients or 
predictors for the statistical weather prediction. This method, called EOF analysis, 
separates the data sets into eigenmodes. Generally speaking, each mode has an associated 
variance, a non-dimensional spatial pattern, and a dimensional time series. From this data 
matrix a 255-square spatial covariance matrix is calculated by 
... 

R,, R^, ... R^^ 

R = 

R^, R^^ ... R^^ 
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where n and ot( 1,2,...,255) denote the grid locations. The diagonal elements of the 
covariance matrix =(n = l,2,...,N) are the variance at location r„ . The off-diagonal 
elements are the covariance with spatial lag equal to the difference between the row and 
column indices. This symmetric matrix has its 255 real eigenvalues , and eigenvectors 
(/>airj) such that 

= i = l2,...N (12) 

M 

The eigenvectors ^,,^2,...,^^ are called Emprical Orthogonal Functions. Each is a 
255-point (17x15 grid) distribution of SST anomaly pattern. The eigenvalues, 
A„(a = 1,2,. ..,Af), are all positive and the summation of them, , equals the total 

variance. Therefore, is considered as the portion of total variance “explained” by the 
EOF (j)^ . It is convenient to label the eigenfunctions so that the eigenvalues are in 
descending order, i.e., 

(13) 

2. Principal EOF Modes 

Intraseasonal and interannual anomalies can be quantitatively investigated with 
the method of emprical orthogonal function (EOF) analysis and a number of its 
generalized forms. In order to delineate the major modes of variability in the winter and 
spring seasons (December to May) of 21 years (1964-84), we perform a conventional 
EOF analysis and obtain the first six leading EOFs. The first six leading EOF modes are 
able to account for 81.7% of the total variance during this season. They are summarized 
in Table 1. When combined, EOFl and EOF2 explain over half (57%) of the total 
variance. Each EOF mode is normalized so that its total spatial variance is equal to unity. 
So, those patterns of the first six EOFs are enough to explain the spatial anomalies of 
whole SCS. Hence, it suffices to focus on the first six EOFs. Their structures are shown 
in Figure 11. Notice that the normalized EOF values are multiplied by 10. 
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EOF 

Variance(%) 

Accumulative Variance 

(%) 

1 

35.6 

35.6 

2 

21.4 

57.0 

3 

12.4 

69.4 

4 

1 __ 1 

74.7 

5 1 

4.7 1 

79.4 


2.3 

81.7 


Table 1. Variances of the First Six Leading EOFs. 


The EOFl mode has an orderly homogeneous pattern of correlation and accounts 
for up to 35.6% of the total spatial variance during the winter and spring seasons 
(December to May). The NE - SW oriented isotherm pattern over this region is believed 
to be related to the ensemble mean pattern. The EOF2 (Figure 1 lb) mode shows a strong 
eddy pattern with a closed isoline of -4°C (actual value is -0.4°C) to the west of Luzon 
Island (116°-119°E, 16°-19°N) and accounts for about 21.4% of variance. The size of 
the eddy is around 100,000 km^. However, we still caimot identify whether the eddy in 
EOF2 represents a warm or cold ring since the values in Figure 8 are normalized. We 
must therefore examine the products of the EOF values and the corresponding time series 
of the amplitudes. 
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3. Temporal Variabilities 

The data matrix, f{r„,T^) is approximately written by 

= ( 14 ) 

where 8^(7^) is a set of time series of size P, representing the temporal variation of the 
associated spatial pattern described by EOF . The result of the time series seems to 

be a projection of the SST anomaly through a ‘filter’ of EOF modes during a time scale. 

The time series of EOFl amplitudes 8,(7^) for 1964-84 are shown in Figure 12. 
The EOFl mode is always negative throughout the whole South China Sea 

(Figure 11a). Therefore, 5,(7^) >0, corresponds to negative SST anomalies, and 
8i(7p)<0, corresponds to positive SST anomalies. Taking the year of 1968 as an 
example, 5,(7^) is always positive from early January to early May, which implies 
negative SST anomalies in the whole SCS during this period. Furthermore, the strength 
of the negative SST anomalies is weak near the west coast of Borneo and increases 
northwestward towards the southeast coast of China (Figure 11 a). 

The time series of EOF2 amplitudes 8^(7^) for 1964-84 are shown in Figure 13. 
The EOF2 mode ^2(6.) (Figure lib) is almost negative throughout the whole SCS 
except near Hainan Island and shows a strong eddy with closed isoline of -4°C to the 
west of Luzon Island. As before, (7^) > 0, corresponds to negative SST anomalies (or a 
cold-core eddy), and B^(7p)<Q), corresponds to positive SST anomalies (or a warm-core 
eddy). Taking the year of 1964 as an example, 8^(7^) is almost always negative from 
early Janutuy to early May and reaches a minimum value (-15) in 1-10 April, which 
implies positive SST anomalies over the whole SCS and a warm-core eddy to the west of 
Luzon Island during the 1964 winter and spring seasons. The strength of the warm-core 
eddy increases from February to April. The SST anomaly of the warm-core eddy in April 
is around 6‘’C (-0.4°C x -15). The actual SST anomaly (6°C) is the product of the 
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second EOF mode (-0.4°C) and its corresponding amplitude (-15). On the other hand, 
taking the year of 1972 as another example, is almost always positive from early 

January to early May with a mean value around 3, which implies negative SST anomalies 
over the whole SCS and a cold-core eddy to the west of Luzon Island during 1972 winter 
and spring seasons. The strength of the cold-core eddy fluctuates from January to May. 
The SST anomaly of the cold-core eddy in 1972 is around -1.2°C (-0.4°C x 3). 

D. SCS COLD-CORE EDDY DURING EL NINO YEARS 

From the above analysis, we see a different thermal behavior in the central SCS 
between the 1964 and 1972 winter monsoon seasons: a warm-core eddy in 1964 and a 
cold-core eddy in 1972. What might cause this interannual variability? Initially, one 
might think of a change in the large-scale surface atmospheric circulation. A reasonable 
guess is the El Nino event (1972 is an El Nino year). During an El Nino year, the 
monsoon circulation was greatly weakened (Kuo, 1989). The wind stress curl over the 
SCS might be different between non-El Nino and El Nino years. Following this line of 
reasoning, we averaged the amplitudes of EOF 1 and EOF2 for El Nino years and non-El 
Nino years, respectively. The set of El Nino years consists of 1965, 1969, 1972, 1976, 
1982, and 1983, and the set of non-El Nino years includes the rest of the years. The mean 
fields for non-El Nino years generally may represent conditions closer to the 
climatological mean. The averaged amplitudes of EOF 1 of non-El Nino years indicate 
positive values from February to March (Figure 14a), which implies negative SST 
anomalies over all the SCS during this period. Furthermore, the strength of the negative 
SST anomalies is weak near the west coast of Borneo and increases northwestward 
towards the southeast coast of China. Conversely, the averaged amplitudes of EOF 1 of El 
Nino years are negative values (Figure 14b). Its SST anomaly pattern is opposite to the 
pattern of the non-EI Nino years during this period of time. 

The averaged amplitudes of EOF2 of non-El Nino years indicate negative values 
during the whole winter and spring seasons (Figure 15a), which implies positive SST 


28 



anomalies over all the SCS and a warm-core eddy to the west of Luzon Island during this 
period. However, the averaged amplitudes of EOF2 during El Nino years indicate 
positive values during the whole winter and spring seasons (Figure 15b), which implies 
negative SST anomalies over all the SCS and a cold-core eddy to the west of Luzon 
Island during this period. The averaged amplitude is usually smaller than individual 
years. This implies that a SCS warm-core eddy appears frequently from January to May 
during non-El Nino years. 

Our data set also shows that during El Nino years a cold-core eddy will frequently 
appear from January to May in the central SCS (Figure 15b). After relating the difference 
in the SST anomaly pattern to non-El Nino and El Nino years, we believe that the 
appearance of a warm-core or a cold-core eddy is related to the global atmospheric 
circulation. 

E. SCENARIOS FOR THE CENTRAL SCS EDDY FORMATION 

The MOODS data set shows that during boreal winter and spring a warm-core 
eddy often appears in the central SCS except for El Nino years when a cold-core eddy 
usually appears. To examine this further, we use climatological data sets such as the 
monthly mean surface wind stress data (Hellerman and Rosenstein, 1983) and monthly 
mean temperature and salinity data sets (Levitus, 1982) to study the eddy formation 
mechanisms in non-El Nino years since the mean features of the non-El Nino years are 
more likely closer to the climatology. 
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Figure 14. Time Series (January to May) of EOFl Amplitudes Averaged Anomaly: 
(a) Non-Nino Years, (b) and El Nino Years (1965,1969, 1972,1976,1982,1983). 
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Figure 15. Time Series (January to May) of EOF2 Amplitudes Averaged Anomaly: 
(a) Non-Nino Years, (b) and El Nino Years (1965,1969,1972, 1976, 1982, 1983). 
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We try to answer the following questions; What is the mechanism for a warm- 
core eddy formation under climatological conditions (more likely representing non-El 
Nino years)? What is the difference of the surface forcing over the SCS between non-El 
Nino and El Nino years? What is the mechanism for a cold-core eddy formation in El 
Nino years? 

Two possible mechanisms may explain the central SCS warm-core eddy 
formation during the winter monsoon season; (1) the effects of downwelling, and (2) the 
Kuroshio intrusion. Figure lOd indicates that the early SCS warming starts at the central 
and south SCS, not at the Luzon Strait. This implies that the Kurushio intrusion is not a 
key factor to cause the central SCS warm-core eddy formation. Therefore, the 
downwelling effect is the predominant mechanism for the central SCSWP formation in 
the boreal winter and spring. 

1. Ekman Pumping in Central SCS 

The monthly climatological wind data (Hellerman and Rosenstein, 1983) over the 
SCS and adjacent areas were used to compute the Ekman pumping velocity (Figures 16- 
19). 


(15) 


Here are the surface wind stresses. Figures 16-19 shows the Ekman pumping 


velocity by the color contours for March-June. The arrows indicate the surface wind 
stress. We see that during March-May (Figure 16-18) a downwelling region with a 
vertical velocity around -0.025 cm/s appears in the central SCS. After the summer 
monsoon onset (June), the surface winds change direction, and a weak upwelling area 
with vertical velocity around 0.005 cm/s appears in the central SCS between 14°N and 
19°N (Figure 19). 
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2. Thermal Feature of the Central SCS 

Figure 20 shows two sets of April MOODS temperature profiles. One is for 1968 
(a non-El Nino year), and the other is for 1982 (an El Nino year). The central SCS water 
is always highly stratified before the summer monsoon onset. 

The ocean mixed layer temperature is much higher than the sublayer temperature. 
The Ekman downwelling in the central SCS during the winter and spring seasons will 
suppress the vertical cold advection from the sublayer. This downwelling effect, plus the 
increase in solar radiation in spring (April-May), may provide some important pre¬ 
conditions for the central SCS warm-core eddy formation. 

3. An Air-Sea Feedback for a Warm-Core Eddy (Non-El Nino) 

We propose an air-sea feedback scenario for the central SCS warm-core eddy 
formation in non-El Nino years (Figure 21). From late winter to spring, an anticyclone 
appears in the central SCS from the ensemble mean atmospheric surface streamline 
analysis (Figure 22a). This anticyclone generates downwelling in the central SCS and in 
turn prevents the cold deep water from being advected to the surface. This will promote 
the formation of the warm-core eddy in the central SCS. As the central SCS warm-core 
eddy persists, being around 1 °C warmer than the surroundings, an atmospheric surface 
low pressure center will be generated above the warm water, which may trigger the 
summer monsoon onset. 

After the summer monsoon onset, an atmospheric surface cyclone will occupy the 
central SCS (Figure 22b). This cyclone will generate upwelling in the central SCS, which 
entrains the deep cold water into the surface mixed layer. This upwelling effect will 
finally destroy the central SCS warm-core eddy. 

During El Nino events, the Asian monsoon circulation is greatly weakened (Kuo, 
1989) and so will be the downwelling effect. Therefore, the warm-core eddy is not likely 
to occur in the boreal winter and spring during El Nino events. The mechanism for a 
cold-core eddy foimation in El Nino years should be further studied. 
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Figure 16. The Ekman Pumping Velocity During March. The Colors Associated 
With Positive (Negative) Values Indicate Upwelling (Downwelling) (cm/s). The 
Arrows Denote the Mean Wind Stresses (From Hellerman and Rosenstein, 1983). 
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Figure 17. The Ekman Pumping Velocity During April. The Colors Associated With 
Positive (Negative) Values Indicate Upwelling (Downwelling) (cm/s). I'he Arrows 
Denote the Mean Wind Stresses (From Hellerman and Rosenstein, 1983). 
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Figure 18. The Ekman Pumping Velocity During May. The Colors Associated With 
Positive (Negative) Values Indicate Upwelling (Downwelling) (cm/s). The Arrows 
Denote the Mean Wind Stresses (From Hellerman and Rosenstein, 1983). 
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Figure 19. The Fkman Pumping Velocity During June. The Colors Associated With 
Positive (Negative) Values Indicate Upwelling (Downwelling) (cm/s). The Arrows 
Denote the Mean Wind Stresses (From Hellerman and Rosenstein, 1983). 
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Figure 21. An Air-Sea Feedback Scenario for the SCS Warm-Core Eddy Life Cycle. 


43 














IV. NUMERICAL SIMULATIONS 


We use a numerical model to simulate the formation of the central SCS warm- 
core eddy. Through this simulation, we hope to verify the downwelling effect scenario. 

A. PRIMITIVE EQUATION MODEL 


The model we use is the three dimensional model developed by Blumberg and 
Mellor (1983, 1987) with the hydrostatic and Boussinesq approximations (Bryan, 1969) 
and has the following features: (1) horizontal curvilinear coordinates and an "Arakawa C" 
scheme (Arakawa and Lamb, 1977), (2) sigma coordinates in the vertical with realistic 
bathymetry, (3) a free surface, (4) a second-order turbulence closure model for the 
vertical viscosity (Mellor and Yamada, 1974, 1982), (5) horizontal difflisivity coefficients 
calculated by the Samagorinsky (1963) parameterization, and (6) split time steps for 
barotropic and baroclinic modes. The basic equations are given by 


(16) 




ae ae d^ 

— + V o V^-i- w— = —(7f„ —) + Fg 

at az az. az 

a az az az ^ 


where V = (u,v) is the horizontal velocity, w the vertical velocity, and g the gravitational 
acceleration. is the characteristic density, p the in situ density, 9 the temperature, S 
the salinity, p the pressure, the vertical eddy diffusivities for turbulent mixing 
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of momentum {K^) and temperature and salinity {K„), {F^,F^,Fg,F^) the horizontal 
mixing terms representing the subgrid scale processes. The density p is computed by the 
equation of state of the form 

P = piG,S) (17) 

given by Fofonoff (1962). 

The model domain is 99°E-121°E, 3°S-25°N which includes the South China 
Sea and the Gulf of Thailand. With its 20 km horizontal resolution and 23 vertical sigma 
coordinate levels, it has 125x162x23 grid points (Figure 23). All depths less than 10 m in 
the domain are set to be 10 m in order to satisfy the constraint H+t]>Q. In the vertical 
the domain is divided into 23 levels cr„(M = 1,2,...,23). Here cr„ = nAcr,Acr = 
the simulation, we use the Samagorinsky parameterization to compute the horizontal 
mixing coefficients with C=0.1, which results in the horizontal viscosity ranging from 
200 to 500 m^s '. The background vertical mixing coefficient is 10^ m^s'' (Blumberg and 
Mellor, 1983, 1987). 

B. LATERAL BOUNDARY CONDITIONS 

There are several major straits connecting the SCS to the Pacific Ocean and 
surrounding seas. Many are very shallow (i.e., the Strait of Malacca) or irregular (i.e., the 
Balabac Strait) and difficult to be adequately represented in a numerical model. 
Therefore, we close the Strait of Malacca, and all the small straits between Luzon Island 
and Borneo in the numerical model. We combine the Karimata Strait and the Gasper 
Strait into one open boundary. After this treatment, our model has three open boundaries: 
a northern boundary (Taiwan Strait), an eastern boundary (Luzon Strait), and a southern 
boundary (Karimata-Gasper Strait). The transports at the three open boundaries were 
determined according to Wyrtki's (1961) estimation with some modification due to the 
boundary treatment. Table 2 shows the seasonal variation of the transport 
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Figure 2: 




Month 

Feb. 

Apr. 

Jun. 

Aug. 

Oct. 

Dec. 

Gasper & Karimata Straits 
(Eastward Positive) 

4.4 

0.0 

-4.0 

-3.0 

1.0 

4.3 

Luzon Strait 
(Eastward Positive) 

-3.5 

0.0 

3.0 

2.5 

-0.6 

-3.4 

Formosa Strait 
(Northwards Positive) 

-0.9 

0.0 

1.0 

0.5 

-0.4 

-0.9 


Table 2. The Bi-monthly Variation of Mass Transport (Sv) at The Open Boundaries. 
The Values were Taken from Wyrtki (1961). 




The barotropic velocity at the open boundaries is obtained by 
bt _ Transport 
" ' \,hdl 

where / is the boundary width, and n means the normal direction to the boundary. The 
baroclinic velocity, K/‘'(c 7) is calculated by 

= (19) 

1 -e 

During the numerical integration, both and F/*" were interpolated into the 
time step from the bi-monthly values of the volume transport listed in Table 2. 

On the lateral open boundaries temperature (T) and salinity (S) are prescribed at 
the inflow boundaries, whereas at the outflow boundaries, the advection equation, 
namely, 

1(7,5)+ F„4(r,S) = 0 (20) 

a ai 

is solved for both barotropical and baroclinic modes. 

C. STABILITY CONSTRAINTS 

The stability constraints were derived from various linearized subsets of equations 
(Blumberg and Mellor, 1987). Since Ax=20 km, Hmax^SOOO m, and Umax = J fnJs, the 
typical CFL condition derived from Blumberg and Mellor (1987) for the external mode is 
45 5 and for the internal mode is 2000 s. In our numerical model, the time step was set as 
25 s for the external mode and as 900 s for the internal mode. 

D. SURFACE FORCING 

The atmospheric forcing includes wind forcing and thermodynamic forcing. We 
took the monthly mean climatological wind stress data (Hellerman and Rosenstein, 1983) 
as the wind data for the middle of the month and interpolated them for each day. The 
wind stress has a typical magnitude of 1 dyne/cm^ (corresponding to a surface wind speed 
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of 10 m/s). Northeasterly winds dominate during the winter and spring seasons while 
southwesterly winds prevail during the summer monsoon season. 

There are two approaches to setting up the surface thermal forcing: (1) flux 
forcing when the surface heat and salt fluxes, m'^ ) and 25'(m s'* ) are given, and 

(2) restoring forcing when the fluxes are not given and the model surface temperature (T) 
and salinity (S) are relaxed to the observed values (TqI,) and (Sob) • 

K^^ = a^Qf,+a^CiT„,-T), (21) 

- S), 

dz 

where /f//is the thermal exchange coefficient and C is the relaxation constant taken to be 
0.7 m/d, which is equivalent to a relaxation time of 43 days for an upper layer 30 m thick, 
(a, ,0(2) are (0,l)-type stvitcher parameters: when a, = \,a-^ = 0, only flux forcing is 
present; a, = 0,0:2 = only th® restoring forcing is present, o:, = 1,0:2 = 1> both flux 
forcing and the restoring forcing are present. Process studies usually use only the 
restoring forcing as the surface thermal forcing (e.g., Brady and Gent, 1994). Therefore, 
in this simulation we use only the restoring forcing. We construct the (7),j ,S„j) data 
through interpolation from the monthly mean temperature and salinity data (Levitus, 
1982). Both the horizontal and vertical resolutions of the Levitus (T,S) data are different 
from our model resolution. Hence, we use optimal interpolation (01) to convert the 
Levitus (T,S) data to our model grids. 

E. MODEL SIMULATED SST FIELDS 

The model was initiated from April mean (T,S) fields (Levitus, 1982) and forced 
by the surface wind stress (Hellerman and Rosenstein, 1983) and thermodynamical fluxes 
depicted in (21). The model was integrated for three years. The third year’s outputs were 
used for discussion. 

Three days of the simulated surface currents were picked for investigating the 
time evolution of the circulation pattern. We choose 30 January for the winter monsoon 
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forcing, 30 March for the boreal spring forcing, and 30 April for the pre-summer 
monsoon forcing. 

On 30 January a strong southward coastal jet is present near the Vietnam coast, 
and no cyclonic or anticyclonic rotation is observed in the central SCS (Figure 24a). The 
isolines of SST in the central SCS are almost parallel to the southeast China coast (Figure 
25a). 

On 30 March the strength of the southward Vietnam coastal current is greatly 
reduced. An anticyclone with a strong northward branch is generated between 13°N- 
18°N, 111 °E-115°E (Figure 24b). Several small cyclones are also simulated, e.g., a small 
and weak one (next to the anticyclone) is located at 13°E-16°E, 115°E-116°E. A warm- 
core eddy is also simulated associate with the anticyclone (Figure 25b). 

On 30 April the southward coastal jet disappears near the Vietnam coast. In the 
SCS an anticyclone persists between 13°N-18°N, 111°E-115°E (Figure 24c). A warm- 
core eddy is also simulated associate with the anticyclone (Figure 25c). 

The simulations confirm our scenario for the warm-core eddy formation under the 
climatological forcing. The cold-core eddy formation during El Nino events should be 
studied in the future research. 
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Figure 24. Model Simulation of Surface Circulation for (a) 30 January, (b) 30 
March, and (c) 30 April. 
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Figure 25. Model Simulation of SST (°C) for (a) 30 January, (b) 30 March, and (c) 
30 April. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

The goal of our study was to detect the meso-scale eddies in the central South 
China Sea during the winter and spring seasons from the Na\'y’s Master Observational 
Oceanographic Data Set (MOODS), and to investigate mechanisms for the eddy 
formation. Optimal interpolation, a composite analysis and an EOF analysis were used to 
study the first objective; a 3-D numerical model was used for the second objective. 

After quality control procedures were employed, the data set was composed of 
180,059 observed profiles spanning the years 1964 to 1984. Due to a sparsity of data, ten- 
day bins were constructed for each year. Optimal interpolation (OI) was used to establish 
a gridded SST data set for the SCS. We applied optimal interpolation to construct a ten- 
day interval synoptic data set (756 total) from 1-10 January 1964 to 21-31 December 
1984 on a 0.5° x 1° grid (finer resolution in the zonal direction) from the MOODS 
temperature data. The gridded data size is 45x25x756. 

The composite analysis on the gridded data set indicates that the annual mean SST 
seasonal variation has four stages; (1) a cooling phase, (2) an early warming phase, (3) a 
central SCS warming phase, and (4) a mature warming phase. A warm-core eddy starts to 
occur in the early warming phase and persists imtil the end of the central SCS warming 
phase. It disappears in the mature warming phase. 

The EOF analysis of the SST anomaly data shows an eddy in the central SCS 
being present in the second principal mode. The eddy has distinctive thermal features 
between non-El Nino and El Nino years. During non-El Nino years the central SCS has a 
warm-core eddy. However, during El Nino years the central SCS has a cold-core eddy. 
This implies that the eddy formation in the central SCS is related to the global 
atmospheric circulation. 


55 




A hypothesis for the SCS eddy formation under the climatological forcing (more 
likely non-El Nino years) was also proposed in this study. During late winter an 
atmospheric anticyclone over the central SCS generates downwelling and in turn prevents 
the cold deep water from being advected to the surface. As solar radiation increased in the 
late winter monsoon,, the central SCS surface water is heated more than the surrounding 
waters by the downwelling effect. This will promote the formation of a warm-core eddy 
in the central SCS. As the central SCS warm-core eddy persists, being a few degrees 
wanner than its surroundings, an atmospheric surface low pressure center will be 
generated above the warm water, which may trigger the summer monsoon onset. After 
the summer monsoon onset, an atmospheric surface cyclone occupies the central SCS. 
This cyclone will generate upwelling in the central SCS, which advects the deep cold 
water into the surface mixed layer. This upwelling effect will finally destroy the central 
SCS warm-core eddy. This downwelling hypothesis on the SCS warm-core eddy 
formation during non-El Nino years was verified by a 3-D primitive equation model. The 
air-sea feedback scenario has not been verified for El Nino events. 

B. RECOMMENDATIONS 

We recommend that further studies should incorporate (1) 3-D thermal features of 
the SCS meso-scale eddies, (2) constructing a composite for the SCS surface wind fields, 
to test if they are upwelling or downwelling favorable wind field; (3) running the 3-D 
numerical model under the El Nino wind forcing. 
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APPENDIX D. TIME SERIES OF EOF 3 AMPLITUDES AVERAGE ANOMALY 


El Nino years average EOF3 



Months(10 days/pt) 
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APPENDIX E. TIME SERIES OF EOF 4 AMPLITUDES AVERAGE AJVOMALY 
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APPENDIX F. TIME SERIES OF EOF 5 AMPLITUDES AVERAGE ANOMALY 
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